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Abstract 

We present quantitative results for the drift velocity of a polymer in a gel if 
a force (e.g. through an electric or magnetic field) acts on a tag, attached to 
one of its ends. This is done by introducing a modification of the Rubinstein- 
Duke model for electrophoresis of DNA. We analyze this modified model with 
exact and Monte Carlo calculations. Tagged magnetophoresis does not show 
band collapse, a phenomenon that limits the applicability of traditional elec- 
trophoresis to short polymers. 
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If electrophoresis is applied for the separation of large polymers, for instance fragments 
of DNA, the velocity becomes independent of the length. This phenomenon is called "band 
collapse" and is widely known amongst experimentalists. A consequence is that long poly- 
mers cannot be separated efficiently as they will all travel with approximately the same 
velocity. Surprisingly, these experimental results can be predicted from a quite simple one- 
dimensional lattice gas model, proposed by Rubinstein [|l| as a model for the dynamics 
of entangled polymers (the repton model), and adapted by Duke as a model for gel elec- 
trophoresis of DNA 1^. It has been studied both analytically and by Monte Carlo 
simulations The same lattice gas model, but seen in a different context, is of inter- 
est in other areas: it describes coarsening , diffusion of hard-core particles [0 with open 
boundary conditions, and interface growth [Q. 

In this work we focus on its application to the motion of long polymers in a gel. We 
study a modification of the traditional electrophoresis: instead of pulling uniformly charged 
polymers with an electric field, we propose to pull on one of the ends of the polymer, which 
can be realized by attaching a magnetic tag to the polymer and applying a magnetic field. 
We will show that according to the repton model, in this case the drift velocity stays length- 
dependent, allowing for the efficient separation of long polymers. In previous work on the 
application of tags in electrophoresis a force acts on the polymer and not on the tag, 
whereas we propose to pull on the tag and not on the polymer. 

First we will give a short description of the repton model, list some of its properties, 
and compare the properties of the repton model with experimental results on DNA elec- 
trophoresis in agarose gels. Next, we will use the repton model to study what happens if a 
force pulls only on one end of the polymer rather than uniformly along the whole molecule. 
Finally, we will outline some possible experimental realizations and summarize our results. 

The repton model is illustrated in Fig. 1. In this model, the gel is represented as a 
collection of pores that are arranged in a square lattice. A polymer, represented as a chain 
of polymer segments or reptons (the black dots in Fig. 1) moves through this gel. 

These reptons move from pore to pore diagonally according to the following rules: 
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(a) reptons in the interior of the chain move only along the sequence of pores occupied by the 
chain. This restriction ensures that the only possible mechanism of movement is reptation. 

(b) at least one repton must remain in each pore along the chain but otherwise the number 
of reptons in a pore is unrestricted. By imposing this rule we achieve that the polymer has 
some elasticity but does not stretch to infinite length. 

(c) the two chain end segments may move to adjacent pores provided that the rule (b) is 
not violated. 

The dynamics of the model are defined as follows: Each possible move (each repton mov- 
ing in any direction) is tried with a constant rate, setting the unit of time. This represents 
the diffusive motion of the polymer. To make it also a model for electrophoresis, Duke has 
introduced a small modification to this original repton model. To each repton a negative 
electric charge is assigned on which the electric field E acts, breaking the symmetry. In- 
stead of unit rates for the allowed moves in the positive and negative x-direction, the rate 
for allowed moves in the positive x-direction (direction opposite to the electric field) equals 
exp(+£'/2), while the rate for allowed moves in the negative x-direction equals exp{—E/2). 

For the application to electrophoresis the quantity of interest is the velocity of the poly- 
mer as a function of the field strength and the length of the polymer. In extensive Monte 
Carlo simulations |^ it was found that a data collapse could be obtained if the rescaled 
velocity vN"^ is plotted as a function of the rescaled electric field NE. Plotted in such a 
way, all data points were found to be located on one single curve, given by 
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where (3 = 5/6 and, as expected from the Nernst-Einstein relation, a/ (3 = DN"^ = 1/3 
where D is the diffusion coefficient of the polymers. For small NE, the rescaled velocity 
vN'^ increases linearly with NE. corresponding to v = E/3N. For large NE, the rescaled 
velocity shows a quadratic behavior. This corresponds to v ^ E'^. 

It turns out that a collapse of experimental data for DNA electrophoresis in agarose can 
be achieved with Eq. 0, as illustrated in Fig. 2 for data points published by Heller et al. |10 



and rearranged as in [|TT|. Taking the velocity v in cm/h, the electric field E in V/cm, and 
DNA lengths in kilobase-pair (kb), one gets a/ (3 = (0.5 ±0.1^/2 and p = (17 ±5)^-^/^ 
where p is the gel concentration in percent agarose by weight. 

This surprisingly good agreement of experimental data with the theoretical prediction (|1]) 
indicates that the repton model contains all important ingredients for DNA electrophoresis: 
entanglement of the polymer in the gel, entropic elasticity of the polymer, and reptation as 
the dominant mechanism of movement. 

In this work we use the repton model to study the behavior of polymers in a gel if a force 
is applied only to one of the ends of the polymers. To do this within the framework of the 
repton model, we preserve the rules that specify which moves are allowed and which ones 
are not. We assume that the force acts only on repton 1. The attempt frequencies are now 
not uniform along the chain: as in the repton model for DNA electrophoresis as adapted by 
Duke, repton 1 moves in the positive and negative x direction with rates q = exp{E/2) and 
= exp{—E/2), respectively; as in the original repton model of Rubinstein, all allowed 
moves of all other reptons occur with unit rate. 

These dynamical rules may be conveniently rephrased in lattice gas language by writing 
a master equation for a one-dimensional exclusion process, where the links of the repton 
chain are mapped to the sites of a one- dimensional lattice with L = A^ — 1 sites. For each 
link i = 1 . . . L (connection between reptons i and z + 1) we distinguish three cases: the 
link has a negative, a zero, or a positive component in the x-direction. We identify these 
three cases with the presence of an A-particle, a vacancy, and a i?-particle, respectively, 
on site i of the lattice. By this mapping the information on the motion of the polymer 
perpendicular to the field gets lost, but this motion is purely diffusive and of no interest for 
the study of the drift velocity of the polymer in direction of the field. The drift velocity is in 
particle language given by the difference between the current of A-particles and the current 
of i?-particles along the chain. 

In terms of particles, the dynamics may be described as follows: If one site of a pair of 
nearest-neighbor sites is occupied by a particle while the other site is empty, this particle 



hops to the empty site with rate 1. In addition to these hopping events, particles can 
be created and annihilated at the ends. At site L, A- and S-particles are created and 
annihilated with rate 1, while at site 1 74-particles are created and annihilated with rates q 
and and S-particles with rates and q, respectively. 

By writing this process in a quantum Hamilton formalism, the stochastic time evolution 
of the system is given by the Hamiltonian H = + 6^(1) +Yli=i 'Wi(l) of a spin-1 quantum 
chain where bi{q) — q{l — nf — af — hi) + ^'"'^(1 — nf — ai — bf), and Ui{q) — q{nfn^^i + 

Heren^ = El\ nf = 

and n° = 1 — — nf = Ef^ are projection operators on states with spin = 1,0, — 1 
resp. on site i, identified with an 74-particle, vacancy and S-particle resp. The operators 
Oj = Ef^, af = El"^, hi = E];^, and hf = Ef^ are annihilation and creation operators for A- 
and B particles. Ej^ is the 3x3 matrix with matrix elements {Ei'^)a^i3 = 8j^a^k,i3 acting on 
site i. 

The average densities {nf'^) at time t satisfy the continuity equations ^{nf'^) — 
— {nf'^H) = {jfli) — {jf'^) with the particle currents 

ijf'^) = {nf'^nl,) - {n^ntf,) l<^<L-l 

{3^) -q{n^) 

Among the quantities of interest are particle densities and correlations in the steady state 
of the system, in particular the (space-independent) stationary particle currents {jf'^) = 
j^'^ which give information on the internal structure of the polymer and which give its 
average drift velocity v — — . Defining — + and rrii — nf — nf , stationarity 
of the probabihty distribution leads to the following relations: 

J = «i) - {<) (1 < i < L - 1) (2) 

= — 2 — ~ )^ — 2 — ^ ' 
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(0) and (^) give the exact average density of vacant sites at site i in terms of J 



= - [1 - (3L + I) J] + Ji 



(8) 



and therefore the probabihty 1 — of finding only one repton in the pore occupied by 
the i*^ repton. C = L — I]f=i("-i') = 2L/3 + J'L{3L — l)/6 gives the average length of the 
polymer in terms of J^. Eqs. (|^), @ and (H) (for i = 1) give 3 linear equations for the four 
quantities v, J', (mi), (n?) relating the currents to the average densities at site 1. This may 
be used to derive an exact expression for v in terms of any of the three other expectation 
values. 

Now we study f as a function of q and N = L + 1. The diffusion coefficient D of the 
repton chains has not changed from the original model, and is to lowest order in given by 
D = 1/ (3A^^) [0. The force on the chain is now independent of its length and proportional 
to E. Combining this with the Nernst-Einstein relation v ^ F ■ D yields v = E/SN"^ for 
small forces that do not disturb the equilibrium configuration of the chain. 

For intermediate fields and short chains, we can obtain the drift velocity numerically 
exactly by computing the eigenvector with eigenvalue of the quantum Hamiltonian H. 
Results for chains with up to 7 reptons {L = 6) are presented in Fig. 3. For intermediate 
fields and longer chains, we study the model by Monte Carlo simulations. We performed 
these simulations by means of a fast "multi-spin" algorithm as described in an earlier pa- 
per 0. Results are included in Fig. 3. 

In the limit of an infinitely strong force, we can solve the model exactly for any L. For 
an infinite creation rate of A-particles and correspondingly infinite annihilation rate of B 
particles on site 1 one has {n^) = 1 and {nf) = (nj) = 0. In this limit the stochastic process 
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becomes equivalent to a process on a chain on L' = L — 1 sites where A particles are created 
(annihilated) with rate 1 (0) on site 2 and where no 5-particles are created or annihilated. 
Going through the same calculations as above gives expressions similar to (0-0), but with 
(D and (D replaced hj J' = v = (ra^) and the range of definition in (0) and (||) replaced 
hj2<i<L~l. Together with (^) (for i = 2) this shows that the velocity approaches 
a constant that depends on the length of the chain: v = 1/(3L — 2) = 1/(3 — 5). This 
asymptotic regime starts around E = 3 a.s can be seen from the exact and Monte Carlo 
results shown in Fig. 3. 

Thus, according to this modified repton model, there is no band collapse if we pull on 
one end of the chain, even in the limit of an infinitely strong field. In this limit, the polymer 
is not fully stretched but has an average length C = L(5 + (3L — 2)^^)/6 ~ 5L/6. 

Experimentally, a realization of pulling only at one end of each polymer, is to attach 
to the end a magnetic bead. Magnetic beads are produced commercially; for a discussion 



of their properties and traditional application in biotechnology we refer to |T3[. These 
commercially available magnetic beads are rather large for use in a normal agarose gel, and 
have the risk to get stuck, but a dense polymer solution may be a suitable replacement for 
the agarose gel. It is interesting to point out that according to our model, variations in 
the force on the end of the polymer, (e.g. due to differences in the magnetic beads) do not 
lead to a significant band widening if the force is sufficiently large, as in that case the drift 
velocity is independent of the force (see figure 3). 

If such "magnetophoresis" turns out to be possible, this experimental approach might 
also avoid band collapse in RNA electrophoresis and make sequencing of longer fragments 
possible. Another interesting topic is the behavior of a polymer that experiences a magnetic 
force on one end, in addition to a uniform electric force on the whole chain. 

Considerable effort has been invested into solving the Rubinstein-Duke model for elec- 
trophoresis, but to the best of our knowledge no rigorous results have been obtained, except 
for very short polymers and for a related model with periodic boundary conditions [Q. 
In zero field the model is equal to the symmetric exclusion process with tagged particles, 



which is (at least) partially integrable by the matrix Bethe ansatz discussed in [jl^. An open 
question is whether open boundary conditions as proposed here leave the system completely 
integrable; this would allow for a full exact solution of both the model proposed here and 
the original Rubinstein model. 

To conclude: previous work has shown that the Rubinstein-Duke model for DNA elec- 
trophoresis provides a good description of experimental results of electrophoresis of DNA 
fragments in agarose gels fill]. We propose a modification in the model to describe pulling 
only on the end of the polymer instead of uniformly along the chain. With Monte Carlo 
simulations and exact calculations we predict that this provides a way of avoiding band 
collapse, a phenomenon that severely limits the separation of long polymers by traditional 
electrophoresis. Experimental realizations of this modified model and various applications 
are discussed. 

We thank J.-L. Viovy, C. Heller and S. Margel for interesting comments on experimental 
realizations. G.M.S. would like to thank C. Godreche for stimulating discussions. This work 
was supported by the EPSRC under Grant No. GR/J78044 (G.T.B.), and under the Human 
Capital and Mobility program of the European Community (G.M.S.). 
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FIGURES 

FIG. 1. Left: Two-dimensional representation of the repton model of a long polymer. The 
reptons (black circles) are located in the pores of the gel (squares of the lattice). Right: projected 
repton model where the reptons are numbered along the chain, and the x-coordinate of repton i is 
plotted as a function of i. The arrows indicate the allowed moves. 

FIG. 2. Experimental data on the rescaled drift velocity vN'^/P as a function of the rescaled 
electric field NE/a. The solid curve is the theoretical prediction by the repton model, see eq. ||. 
The data points are obtained from and rearranged as in |lT[| . 

FIG. 3. Drift velocity ?; as a function of the force E. The solid lines are exact results for 
L = 2 . . . 6, the dashed lines are obtained from Monte Carlo simulations for L=7,10,15,20,30 and 
50. Statistical errors smaller than 2% are omitted. 
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